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V^ , Abstract 

j-^ , Recently reported infrared (IR) galaxy number counts and cosmic infrared background (CIRB) all 

(^ • suggest that galaxies have experienced a strong evolution sometime in their lifetime. We statistically 

f^ , estimate the galaxy evolution history from these data. The evolution of galaxy luminosity is treated as a 

^D ' nonparametric stepwise function, in order to explore the most suitable evolution history which satisfies the 

r"| . constraint from CIRB. We find that an order of magnitude increase of the far-infrared (FIR) luminosity at 

Q^' redshift z = 0.5-1.0 is necessary to reproduce the very high CIRB intensity at 140 /xm reported by Hauser 

JL. ! et al. (1998). We note that an evolutionary factor larger than ten at high-z {z > 2) overpredicts the CIRB 

}^ • intensity at submillimeter wavelength regime. Basically the evolutionary factor increases by a factor of 30 

C/3 , up to z ~ 0.75 and decreases to, even at most, a factor of 10 toward z ^ 5, though many variants are 

Cu ' allowed within these constraints. This evolution history also satisfies the constraints from the galaxy number 

^ i counts obtained by IRAS, ISO and, roughly, SCUBA. The rapid evolution of the comoving IR luminosity 

density required from the CIRB well reproduces the very steep slope of galaxy number counts obtained 

by ISO. This result is also consistent with some recent results of multiwavelength surveys and follow-up 

C^ ■ observations of ISO extragalactic sources. We also estimate the cosmic star formation history (SFH) from 

the obtained FIR luminosity density, considering the effect of the metal enrichment in galaxies. The derived 

SFH increases steeply with redshift in 0<z<0.75, and becomes flat or even declines at z > 0.75. This is 

consistent with the SFH estimated from the reported ultraviolet luminosity density. In addition, we present 

the performance of the Japanese ASTRO-F FIR galaxy survey. We show the expected number counts in 

the survey. We also evaluate how large a sky area is necessary to derive a secure information of galaxy 

evolution up to z ~ 1 from the survey, and find that at least 50 - 300 deg^ is required. 

Key words: Galaxies: active — Galaxies: evolution — Galaxies: formation — Galaxies: starburst — 
Infrared: galaxies 

1. INTRODUCTION frared (hereafter IR) and submillimeter (sub-mm) wave- 

band observations of galaxies have a crucial importance 
The galaxy evolution has long been a strong driving for full understanding of their evolutionary status. 
force of the cosmological studies, and many problems ^^^^^^ -^^^^^^^ ^^^ ^^^^^_^ ^^^^^^^^ ^^^^^1^^ ^^^^ 

have been still unsolved. Not only optical but also in- ^^^^^ ^^^^^^ ^^ ^^1^^^ ^^^^^^ ^^^^^^ compared with that 

Z Research Fellows of the Japan Society for the Promotion of expected from the no-evolution model, and provided a 

Science. new impetus to the related field (e.g. Kawara et al. 1998; 
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frared satellite which will be launched in 2004. De- 
tailed informations of the ASTRO-F mission are avail- 



ablc in http://www.ir.astro.isas.ac.jp/ASTRO-F/index- 



e.html , In order to calculate the expected detection num- 
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Fig. 1. — The assumed galaxy spectral energy dis- 
tribution in the near infrared to radio wavelengths, 
constructed from the color-luminosity evolution in 
IRAS galaxies and the tight correlation of the far-infrared 
(FIR) and radio fluxes. The prominent emission bands 
are PAH features. Vertical dotted lines represent the 
wavelengths of IRAS four bandpasses. The SEDs with 
the FIR luminosity of 10^ 10^ 10^^ 10", 10^^^ lO^^^ 
and 10^^ L(T, are shown from the bottom in this order. 



Puget et al. 1999; Dole et al. 2000; Oliver et al. 2000a; 
Kawara et al. 2000; Okuda et al. 2000). The ITO-^im 
slope proved to be d log TV/ d log 5* < — 2.5 at 5* ~ 0.5 
Jy in these new deep surveys, while the slope would be 
~ —1.5 in the no evolution case. Such an excess of galaxy 
number count is generally understood as a consequence 
of a strong galaxy evolution, i.e. a rapid change of the 
star formation rate in galaxies. 

Another important related issue is the cosmic infrared 
background (CIRB), which is the integrated IR light from 
galaxies, especially from those are so faint that are un- 
able to be resolved. Therefore the CIRB provides an 
important information on the past star formation his- 
tory of galaxies including the inaccessible sources. The 
CIRB has detected by COBE (Puget et al. 1996; Fixsen 
et al. 1998; Hauser et al. 1998). The reported CIRB in- 
tensity was surprisingly high {vL^ = 25 ± 7 nWm^^sr^^ 
at 140 ^m and i^L„ = 14 ± 3 nWm^^sr^^ at 240 ^m: 
Hauser et al. 1998), and provided a profound problem 
awaiting to be solved. 

Now next generation space infrared facilities (e.g. 
ASTRO-F, SmTF, NGST and FIRST) are sched- 
uled, and we can expect a vast amount of novel 
observational knowledge with an unprecedented pre- 
cision. Among the above facilities, ASTRO-F (In- 
frared Imaging Surveyor: IRIS) is the Japanese in- 



ber in the survey planned in the ASTRO-F project, we 
have made a simple empirical model for the galaxy num- 
ber count (Takeuchi et al. 1999; Hirashita et al. 1999). 
The applied model was based on the IRAS surveys. 
Therefore, though our previous model has successfully 
worked, it cannot reproduce the observed slopes of the 
novel ISO source counts, and the evolutionary effect is 
no longer satisfactorily estimated. Thus, we need an im- 
proved model to study the detailed observational plans 
and follow-up strategies in other wavelengths. 

There are a number of attempts to predict the IR- 
sub-mm source counts or CIRB properties. Modeling 
methods of the source counts and the CIRB roughly 
fall into two categories: One is the backward approach, 
which is based on the local far-IR (FIR) luminosity func- 
tion (LF) and the observed spectral energy distribution 
(SED) of galaxies in the FIR - sub-mm, with the as- 
sumptions of simple functional forms for the evolution 
(e.g. Beichman, Helou 1991; Pearson, Rowan-Robinson 
1996; Malkan, Stecker 1998; Xu et al. 1998; Takeuchi et 
al 1999). The other is the forward approach, which is 
based on the models constructed by detailed processes 
related to the evolution of galaxies with a number of 
parameters (e.g. Franceschini et al. 1994; Guiderdoni et 
al. 1998). Tan, Silk, Balland (1999) presented an inter- 
esting approach to combine the two methods. These two 
different approaches have their own merits and demerits, 
and thus they must interplay with each other. 

In this paper, we apply the backward empirical ap- 
proach to derive the evolution of FIR properties of galax- 
ies, but assume no specific functional form for the evo- 
lution history. Instead, we treat the galaxy evolution as 
a stepwise nonparametric form and search for the most 
plausible solutions, which well reproduce both of the the 
available IR galaxy number counts and the remarkable 
spectrum of the CIRB in a consistent manner. 

The rest of this paper is organized as follows: We give 
the model description and formulation in Section 2. We 
illustrate the method of our analysis in Section 3. In 
Section 4 we show the most probable solutions to ex- 
plain the observed IR number counts and CIRB. Some 
discussions are presented in Section 5. Section 6 is de- 
voted to our conclusions. We present the expected num- 
ber counts of the ASTRO-F FIR all-sky survey in the 
Appendix. Throughout this paper we use the follow- 
ing cosmological parameter set, unless otherwise stated: 
Ho = 75kms"iMpc"\ go = 0.1, and Aq = 0. 
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Fig. 2. — Upper panel: the applied 60-/im luminosity 
function (LF) in this paper. This LF is derived by Soifer, 
Houck, Neugebauer (1987). Lower panel: a schematic 
representation of its evolution as a function of redshift. 

2. Model Description 

The galaxy number count model is represented by 
SED, LF, cosmology, and galaxy evolution. We review 
the framework of the number count model and discuss 
these ingredients for our model in this section. 

The infrared - radio galaxy SED is modeled by the 
superposition of the following components. For the in- 
frared - sub- mm component, we consider PAH (poly- 
cyclic aromatic hydrocarbon), graphite and silicate dust 
spectra (Dwek et al. 1997). In the case of the radio-quiet 
sources, the emission at millimeter wavelength regime is 
dominated by the synchrotron radiation explained by su- 
pernova remnants (Condon 1992). 

By considering the above, we start to construct the in- 
frared model SEDs based on the IRAS color-luminosity 
relation (Smith et al. 1987; Soifer, Neugebauer 1991). 
Smith et al. (1987) and Soifer, Neugebauer (1991) found 
that the IRAS FIR colors and FIR luminosity are tightly 
correlated. The FIR color-60-/im luminosity relation de- 
rived by Smith et al. (1987) is as follows: 



[/im], and Lqq [Lq] is the intrinsic luminosity evaluated 
at the 60-/im bandpass, Leo = i^L^, at 60 fim. The re- 
lation reported by Soifer, Neugebauer (1991) is shghtly 
nonlinear, but also a monotonic function of Leo- We in- 
terpreted this relation to the dust temperature Tdust~L6o 
relation and calculated the modified blackbody contin- 
uum with the corresponding Tdust- We adopted the dust 
emissivity Si, oc v"^ with 7 = 1.0, 1.5, and 2.0, and found 
the three values yield only a small difference ( < 10%) in 
the Tdust^Leo relation. We hereinafter use 7 = 1.5. 

We then utilized the data provided by Soifer, Neuge- 
bauer (1991) and derived the approximately linear rela- 
tion between IRAS flux densities, 5*60 and S'25. The tight 
correlation between the FIR and mid-IR (MIR) flux is 
also found in the diffuse emission from the Galaxy by 
Infrared Telescope in Space {IRTS: Shibai, Okumura, 
Onaka 2000). We thus added the mid-infrared spec- 
tra proposed by Dwek et al. (1997) to the FIR com- 
ponent so that the superposed spectra reproduce the 
correlation reported by Smith et al. (1987), such that 
log 525/5*60 — —0.9. The unidentified infrared bands 
(UIBs), which we assumed to be produced by PAHs, are 
also an important component of the IR SEDs of galaxies. 
We set the continuum-to-band intensity ratio as reported 
by Dwek et al. (1997). Detailed properties of PAHs are 
taken from AUamandola et al. (1989), e.g. we set the 
PAH features at 3.3 /im, 6.2 /im, 7.7 /im with a broader 
component at 5.5 — 9.5 /im, and 11.3 /im. 

A remarkably tight and ubiquitous correlation is well- 
known between the FIR continuum flux and radio con- 
tinuum flux (e.g. Helou, Soifer, Rowan-Robinson 1985; 
Bregman, Hogg, Roberts 1992; Condon 1992). For the 
longer wavelength regime, a power-law continuum pro- 
duced by synchrotron radiation (oc 1'^°') dominates the 
observed emission. We set a = 0.7 according to Condon 
(1992), and added the FIR composite spectra, using the 
relation 
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(Condon 1992). Equation (||) is adopted in Helou et 
al. (1988) and will be apphed again in Section 5. This 
is the final SED we use in our number count and CIRB 
models (see Figure 1). 

2.1. Local Luminosity Function and Evolutionary Effect 
We adopted the 60-/im LF based on the IRAS by 
Soifer et al. (1987) as the local IR LF of galaxies: 



where S\ is the detected flux density at wavelength A log (/)o (-^60 [Lq 
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where (Jjq is the number density of galaxies in 
Mpc~'^dex~^. We show the LF in the upper panel of Fig- 
ure 2. We applied the double power-law form for the local 
LF. The lower cut-off luminosity do not seriously affect 
the result, as stated in Takeuchi et al. (1999). The effect 
of the faint-end slope of the LF on the number count es- 
timation is briefly summarized in Takeuchi, Shibai, Ishii 
(2000a). 

We assumed a pure luminosity evolution in this study. 
In this case 60-/Ltm luminosity of a certain galaxy at red- 
shift z is described as 



Leo{z)^Leo{0)f{z) 



(4) 



We also assumed that the luminosity evolution is 'uni- 
versal', i.e. independent of galaxy luminosity. This is 
depicted in the lower panel of Figure 2. In the previous 
backward approach, a certain functional form f{z) was 
often assumed for the evolution of galaxies, and model 
counts are calculated and compared with observations. 
We, on the contrary, treat the evolution of galaxy lumi- 
nosities as a stepwise nonparametric form which is to be 
statistically estimated, in order to explore the most suit- 
able evolutionary history which reproduces the present 
observational results. Mathematically, we regard the evo- 
lution estimation from the number count as an inverse 
problem, and we apply the universal LF evolution in or- 
der to obtain a statistically stable solution. Constraints 
of f{z) from the observed CIRB spectrum will be given in 
Section 4.1. Another approach mainly based on the den- 
sity evolution is discussed in Pearson et al. (2000). We 
note that the present quality of the IR number count data 
is not sufficient to estimate the evolution of the shape of 
the LF, and fine tuning of the LF shape is less meaning- 
ful. Such works remain to be solved after the forthcoming 
release of huge databases. 

2. 2. Formulation of the Galaxy Number Count 

Using the above formulae, now we calculate the flux- 
number (log TV- log S") relation, or the so-called galaxy 
number count. We assume that galaxies can be re- 
garded as point sources (i.e. the cosmological dimming of 
the galaxy surface brightness is not taken into account). 
Then the relation between the observed flux S{h') and the 
emitted monochromatic luminosity L(t'cm) = L{{1 + z)]^) 
is given by 
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where c^l is the luminosity distance. When we fix a cer- 
tain S{i'), we obtain L((l + z)^) by using equation (ra). 
Then the corresponding L6o(S'(^), z) at a redshift z is 
uniquely determined. We define N{> S{h')) as the num- 
ber of galaxies with a detected flux density larger than 
S{v) in a survey solid angle Q, then it is formulated as 



N{> Sii^)) = 



an 

n Jo 



dz 
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dzdnJL,o(s(,. 



{z, L'eo) dL'eo , (6) 
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where d^V /dzdfl is the comoving volume element per sr, 
and Zmax is the maximum redshift we consider in this 
study. As we see later, we set ^max = 5. 

We then formulate the number count with an evolu- 
tion. If a pure luminosity evolution takes place, by us- 
ing equation (0), the evolution of the luminosity function 
with redshift is expressed as 



(z, Leo) dLeo = 0o ( j^j d (^-^ 



(7) 



where 0o(-^6o) is the IR LF in the local Universe. The 
expected number count is expressed as 
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where Leo = L^o/ .f{z). The full formulation of the galaxy 
number count including a density evolution is given in 
e.g. Gardner (1998) or Takeuchi et al. (1999). In this 
paper we call f{z) 'the evolutionary factor'. 

2. 3. Formulation of the Cosm,ic Infrared Background 

The CIRB is generated from the integrated light of 
galaxies. Therefore, combining the SEDs of galaxies and 
number count predictions, we obtain the CIRB spectrum. 
The observed flux density of a galaxy whose IR luminos- 
ity is Leo, S{iy, Lgo), is given by equation (m as follows: 



Siv,Leo) = 



{l + z)L{iy{l + z),L6o) 



indl 



(9) 



where L{iy, Leo) is the monochromatic luminosity of a 
galaxy with Leo. Then the CIRB spectrum /(j^), i.e., 
the background flux density from a unit solid angle, is 
expressed as 



IM 



dz 
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We can deal with the evolutionary effect on I{v) through 
(/)(z,Leo), just the same as in the case of the number 
count calculation. 
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3. Analysis 

In principle the galaxy number count and the spec- 
tral shape of the CIRB are integrated values along with 
redshift. Therefore, the information of galaxy evolution 
with redshift is degenerate if we analyze a dataset ob- 
tained from a single wavelength. But the redshift degen- 
eracy can be solved by treating the multiband observa- 
tional results at the same time as suggested by Takeuchi 
et al. (2000b) , who studied the changes of the IR number 
counts against some extreme galaxy evolution histories. 

Here we thoroughly survey the response of the multi- 
band number counts and the CIRB to the change of the 
evolutionary factor at various redshifts. First we divide 
the redshift into the foUwoing eight intervals, [0 - 0.25], 
[0.25 - 0.5], [0.5 - 0.75], [0.75 - 1.0], [1.0 - 1.5], [1.5 ~ 
2.0], [2.0 - 3.0], and [3.0 - 5.0]. If we set the evolutionary 
factor f{z) = 10 at each interval, the number counts and 
the CIRB will deviate from the no-evolution prediction 
differently, corresponding to each interval. This different 
response of the predictions provides a useful information 
of the evolution with redshift. For example, consider a 
certain deviation in the number count at one waveband. 
Without any other informations, we cannot distinguish a 
certain evolution history from all the other possibilities. 
But when we use the deviations in other wavelengths and 
in the CIRB, we can check the validity of an evolution by 
the goodness of fit to the data at other wavelengths and 
may infer a much more constrained parameter range. 

We show the deviations of the CIRB spectrum and 
the galaxy number counts with an order-of-magnitude 
evolution at each redshift interval in Figures 3 and 4. 
The left panel of Figure 3 depicts the response of the 
CIRB spectrum for hypothetical, an order-of-magnitude 
evolution at the above redshift interval at z < 1. The 
no-evolution prediction of the CIRB is presented by the 
thick solid curve. The dotted line represents the CIRB 
when f{z) = 10 at < z < 0.25, the dashed line at 0.25 < 
z < 0.5, the dot-dashed line at 0.5 < z < 0.75, and the 
dot-dot-dot-dashed line at 0.75 < z < 1.0, respectively, in 
Figure 3. The peak of each CIRB model prediction shifts 
toward longer wavelengths, with the increasing redshift of 
the input test evolution, while the intensity of the CIRB 
at sub- mm regime is insensitive to the evolution at z < 1. 
On the other hand, the right panel of Figure 3 illustrates 
the response of the CIRB caused by each test evolution 
at z > 1. This time the dotted line represents the CIRB 
when f(z) = 10 at 1.0 < z < 1.5, the dashed line at 1.5 < 
z < 2.0, the dot-dashed line at 2.0 < z < 3.0, and the 
dot-dot-dot-dashed line at 3.0 < z < 5.0, respectively. 
Again the thick solid line is the no-evolution prediction. 
At z > 1, the evolution mainly affects the intensity of 
the sub-mni CIRB spectrum. 

How about the number count? We show the effect of 
the test evolution to the multiband number counts in 



Figure 4. The response galaxy number counts at 15 /im, 
60 /im, and 90 /xm are displayed in Figure 4a, while those 
at 170 /im, 450 /x, and 850 /xm are in Figure 4b. The same 
as in Figure 3, the response of the number counts for the 
input test evolution at z < 1 are indicated in the left 
panels, and that for the input evolutions at z > 1 in the 
right panels, in Figure 4a and 4b. 

We clearly see that the MIR 15-/im count is almost 
independent of the status of galaxies at z>0.75, and it 
is quite sensitive to the evolution at 0.25 < z < 0.5. 
Thus the 15-/xm galaxy number count is a useful and 
unique indicator of the evolution at such low redshifts. 
The 60-/xm count depends still strongly on the evolution 
at 0.25 < z < 0.5, and is more sensitive to the evolu- 
tionary status at 0.5 < z < 0.75. On the contrary it 
is insensitive to that at z>l. This trend of the red- 
shift dependence becomes more prominent at FIR 90- 
and 170-/im number counts. The effect of the evolution 
at 1.0 < z < 1.5 begins to appear in the faintest slope 
of the number counts at these wavelengths. The sub-mm 
450- and 850-/xm counts shows, in contrast, a weaker de- 
pendence on the low-z (z < 1) galaxy evolution. They 
are very strongly dependent on the evolutionary factor 
at z > 1. Here, notice that the effect of the test evolution 
at z > 1 are almost constant in 850 /im. This is caused 
by the counter-intuitive, negative X-correction. 

In summary, we are able to resolve the degeneracy of 
the galaxy evolution history by treating the multiband 
datasets and the CIRB at the same time. Since the am- 
plitude of the evolutionary factor at various redshift in- 
tervals correlates with each other, this statistical infer- 
ence problem is still more difficult than a mere fitting 
problem. But we can obtain a set of permitted ranges of 
the evolutionary factors as a function of redshift by this 
method. 



4. Results 

In Figure 5, we summarize the observational con- 
straints of the evolutionary factor estimated from the IR 
galaxy number counts and the CIRB. Thick dashed hor- 
izontal lines depict the bounds mainly based on the IR 
galaxy number counts. Thick solid lines are the con- 
straints based on both CIRB and IR number counts. 
Thick dot-dashed line is determined by the severe con- 
straints from CIRB spectrum at sub-mm wavelengths. 
We set that there are no galaxies at z > 5. This is not a 
strict observational constraints, but the observed CIRB 
suggests that there are few luminous IR galaxies at such 
a high redshift. We explain and discuss the details of 
these constraints in the following. 
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Fig. 3. — The deviations of the CIRB spectrum and the galaxy number counts with an order-of-magnitude evolution 
at the following eight intervals: [0 - 0.25], [0.25 - 0.5], [0.5 - 0.75], [0.75 - 1.0], [1.0 - 1.5], [1.5 - 2.0], [2.0 - 3.0], 
and [3.0 - 5.0]. The left panel depicts the response of the CIRB spectrum for hypothetical, an order-of-magnitude 
evolution at the above redshift interval at 2; < 1. The no-evolution prediction of the CIRB is presented by the thick 
solid curve. The right panel illustrates the response caused by each test evolution at z > 1. 



4-.1. Constraint from Cosmic Infrared Background 

Today we have a considerable information on the spec- 
trum of the CIRB. As we mentioned in Section 1, the 
CIRB has been revealed to have a surprisingly high inten- 
sity at FIR wavelength regime {vL^ — 25±7 nWm~^sr~^ 
at 140 /xm and vL^, = 14 ± 3 nWm^^sr"^ at 240 ^m: 
Hauser et al. 1998). Thus the spectral shape of the CIRB 
can be used to constrain the evolutionary history of the 
cosmic IR luminosity density which is closely related to 
the cosmic history of star formation and metal produc- 
tion. 

We summarize the reported intensity of the CIRB in 
Figure 6. Open squares with downward arrows rep- 
resent the "dark sky" upper limit to the CIRB mea- 
sured by COBE Diffuse Infrared Background Experiment 
(DIRBE). The DIRBE sky brightness varies roughly si- 
nusoidally over the year, due to the complex features of 
the interplanetary dust cloud. The "dark sky" bright- 
ness is that of the darkest area on the sky at each wave- 
length. Open squares with upward arrows represent the 
residual signal of DIRBE after removing the contribu- 
tions from the model foreground sources such as Galac- 
tic diffuse emission, and interplanetary dust emission. 
Filled squares with error bars are the DIRBE detec- 
tions at 140 /im and 240 yum by Hauser et al. (1998). 
The hatched region is the permissible range derived from 
the COBE Far Infrared Absolute Spectrophotometer (FI- 
RAS) high-frequency data, after removal of the cosmic 



microwave background (CMB) signal (Fixsen et al. 1998). 
At wavelength shorter than 140 ^m, the strong emission 
from the interplanetary dust prevents us from estimat- 
ing the reliable value of the extragalactic component, and 
therefore in this study, we do not use the intensity at this 
wavelengths recently reported by some authors. The de- 
tector noise level of DIRBE becomes worse at wavelength 
longer than 240 iim.. On the other hand, the CIRB spec- 
trum at sub-mm regime detected by FIRAS is considered 
to be highly robust, because in this wavelength regime, 
what to be subtracted is only the blackbody spectrum of 
the CMB, which has no model ambiguity. 

Lagache et al. (1999) claim that, from their observa- 
tion of the sky in and around the Lockman Hole, there 
exists a significant contribution from the Galactic warm 
interstellar medium (WIM) to the IR background inten- 
sity. Thus if we adopt their interpretation, the CIRB 
intensities become 15.3 ± 6.4 nWm~^sr~'^ (140 ^m) and 
11.4 ± 1.9 nWm^^sr"^ (240 ^m). The contribution of 
the WIM at sub-mm wavelength is small, and the inten- 
sity of the CIRB at 200 — 2000 /xm is almost unchanged 
from Fixsen et al.'s value. However, in other regions, 
Lagache et al. (2000) reported similar results to Hauser 
et al. (1998) even though they again consider the contri- 
butions of the WIM, and the existence of the WIM still 
seems a matter of debate. We present the claimed value 
of Lagache et al. (1999) by filled triangles with errors in 
Figure 6. 
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Fig. 4a. — The effect of the test evolution to the multiband number counts. The response galaxy number counts at 
15 fim, 60 fim, and 90 /zm are displayed. The same as in Figure 3, the response of the number counts for the input 
test evolution at z < 1 are indicated in the left panels, and that for the input evolutions at z > 1 in the right panels. 



At a glance of Figure 6, it is obvious that no evolution 
prediction is strongly ruled out. The CIRB intensity at 
140 /im suggests an order-of-magnitude evolution of IR 
galaxies. It is important that such a short wavelength of 
the CIRB peak can only be reproduced by a rapid evo- 



lution of galaxies at low redshift {z < 1). Even if we 
use the slightly lower value of Lagache et al. (1999), this 
conclusion is not affected. If we try to explain the peak 
intensity by high-redshift IR galaxies, vast numbers of 
galaxies is required. Furthermore, if in case, the CIRB 
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Fig. 4b. — The same as Figure 4a, except that these counts are at 170 ^m, 450 fi, and 850 ^ni. 



intensity at longer wavelength have to be overpredicted, 
and the observed spectrum will be seriously violated. As 
we noted that the FIRAS spectrum at sub-mm wave- 
length is highly reliable, all of the model must reproduce 
it exactly. At least, models which overestimate FIRAS 
spectrum should be dismissed. The sub-mm slope of the 
CIRB is shallower compared with a single galaxy SED, 



which should be interpreted as an integrated SED of a 
significant numbers of galaxies at high redshift (z > 2 — 3) 
(sec Figure 3). Thus, though the constraint is strong, it 
is not inconsistent with the actual existence of obscured 
high- 2; starbursts (e.g. Hughes 2000; Ishii et al. 2000). 
What we should stress is that there is an observational 
upper bound of the evolution of IR galaxies. In Fig- 
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Fig. 5. — The summary of the constraints of the galaxy 
luminosity evolution inferred from the infrared galaxy 
number counts and the cosmic infrared background radi- 
ation spectrum. 



ure 5, the constraints from the CIRB peak intensity is 
presented as the enhance of the evolutionary factor at 
0.5 < z < 0.75, and the constraints from FIRAS obser- 
vation is shown by the upper bound at 1 < z < 5. 

Hereafter, within the permitted range derived from the 
CIRB, we focus on the three representative evolutionary 
histories. We show these evolutions in Figure 7. Evo- 
lution 1 has a steep rise at < z < 0.5, reaches the 
summit at which the evolutionary factor f{z) is 10, and 
has a plateau at 0.5 < z < 2. Then, at 2 < z, the evolu- 
tionary factor slowly decreases but still higher than the 
local value. This rather moderate evolution is consistent 
with the CIRB 140-/im intensity when we consider the 
WIM contribution, and slightly overestimates the FIRAS 
sub-mm spectrum. These features are also described in 
Figure 6 as the dot-dashed line. 

Evolution 2 has a stronger rise of the evolutionary fac- 
tor than Evolution 1 at < z < 0.5. Beyond the peak at 
0.5 < z < 0.75, the evolutionary factor decreases rapidly 
to the local value at z < 2, and at 2 < z < 5 the IR 
LF is the same as the local one. This evolution mitigates 
the overestimation of the sub-mm CIRB, and the 140-/im 
peak intensity has an intermediate value between Hauser 
et al. (1998) and Lagache et al. (1999). 

Finally Evolution 3 has the strongest rise at < z < 
0.5. This evolution model has a prominent peak at 0.5 < 
z < 0.75, which well reproduces the 140-/im peak of the 
CIRB by Hauser et al. (1998). This evolution also has a 
long plateau at 0.75 < z < 5. The predicted CIRB by 
this model is quite consistent with the FIRAS spectrum. 
Since we present the evolutionary history as a function 
of redshift, we would have an impression that the rise of 
the evolutionary factor at 0.5 < z < 0.75 is unnaturally 
steep, but the peak height is 30 and the duration is ^ 2 
Gyr, both of which are within the bounds of possibility 
and worth being considered. 

These CIRB model behavior is also shown in Figure 6. 
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Fig. 6. — The cosmic infrared background observations 
and our model predictions. Dotted line describes the no- 
evolution prediction. Dot-dashed line shows the model 
CIRB with Evolution 1 , dashed line represents the CIRB 
with Evolution 2, and solid line, with Evolution 3, re- 
spectively. 



No-evolution prediction is represented by dotted line in 
Figure 6. The dot-dashed line is the model CIRB pro- 
duced by Evolution 1, the dashed line is produced by 
Evolution 2, and the solid line, by Evolution 3. 



J^.2. Constraint from IR Galaxy Number Counts 

We examine whether the three evolution models are 
consistent with observed galaxy number counts. Re- 
cently, new observational results of the galaxy number 
counts have become available at MIR, FIR, and sub- 
mm, mainly by the success of ISO (Infrared Space Ob- 
servatory) and SCUBA (Submillimcter Common-User 
Bolometer Array for the James Clerk Maxwell Tele- 
scope). We are able to compile these data as well as 
previously obtained IRAS databases. 

We present the galaxy number counts at 15 /xm - 
850 /im in Figure 8. Again the no-evolution predictions 
are represented by dotted lines, the dot-dashed lines are 
the model number counts produced by Evolution 1, the 
dashed lines are produced by Evolution 2, and the solid 
line by Evolution 3, respectively. Here we go into the 
details of the galaxy count at each wavelength in IR - 
sub-mm. 
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Fig. 7. — Three representative galaxy evolutionary histo- 
ries permitted by the observational constraints presented 
in Figure 5. 



4-2.1. 15-^m galaxy number counts 

The top-left panel of Figure 8 shows the 15-/im galaxy 
counts predicted by our model with Evolutions 1,2, and 
3. We also show the available observed 15-/im ISO- 
CAM galaxy counts (Flores et al. 1999a, b; Clements 
et al. 1999; Aussel et al. 1999; Ahieri et al. 1999; El- 
baz et al. 1999; Oliver et al. 2000a) and the constraint 
from fluctuation analysis {P{D) analysis) of the Hub- 
ble Deep Field (HDF) (Oliver et al. 1997), together 
with IRAS galaxy count (Rush, Malkan, Spinoglio 1993). 
The observed counts are well reproduced by any of our 
three evolutionary histories, and clearly different from 
the no-evolution prediction. The slope and normal- 
ization of the IRAS count must be reproduced by the 
model without any corrections. The ELAIS (European 
Large-Area ISO Survey: Oliver et al. 2000b; for the 



present status, see http://athena.ph.ic.ac.uk) count (Ser- 
jeant et al. 2000) severely constrain the range of the evo- 
lutionary factor, and this shows that the rapid rise at 
z < 0.75 is necessary to explain the MIR galaxy number 
counts. Of course too large evolutionary factor is also 
forbidden by this data. The permitted range of the evo- 
lutionary factor is f{z — 0.5-0.75) = 30 ± 3 if we take 
both the Poisson fluctuation and the cosmic variance into 
account (see Appendix). 



4.2.2. Far-infrared galaxy number counts 

Next we examine the galaxy counts at FIR wave- 
lengths. The top-right panel of Figure 8 is the IRAS 60- 
fj,m galaxy number count (Rowan- Robinson et al. 1991). 
Our model calculations are based on the statistical stud- 
ies on the local IRAS galaxies and the IRAS 60-/im 
LF, the 60-/im count is automatically reproduced at the 
brightest end. On the other hand, the count at the 
faint flux limit restricts the evolution at lowest redshifts 
(z ^ 0.1-0.3). Some previous authors claimed the ev- 
idence of the FIR galaxy evolution at the flux limit 
of IRAS (e.g. Saunders et al. 1990; Ashby et al. 1996; 
Berlin, Dennefeld, Moshir 1997; Springel, White 1998), 
but the reported strengths of the evolution were not con- 
sistent with each other. We found that the evolutionary 
factor can be at most f{z) < 3 in the redshift range 
0.25 < z < 0.5, otherwise the model number count vio- 
late the 1-cr variance range of the observed IRAS count. 

At 90 and 170 /im, deep galaxy counts are obtained by 
ISOPHOT. The photometric calibration of the ISO is stiU 
a matter of debate, and we should take care of the inter- 
pretation of the counts at these wavelengths. The middle- 
left panel presents our model 90-/im galaxy counts with 
observed FIRBACK source counts (Dole et al. 2000) and 
ISOPHOT CIRB project (Juvela, Mattila, Lemke 2000), 
and the constraint from the fluctuation analysis of the 
ISO Lockman Hole data (Matsuhara et al. 2000). We also 
show the IRAS 100 /im galaxy counts by a hatched thin 
area at brightest flux. Again our model shows an excel- 
lent agreement with IRAS count. The middle-right panel 
is the 170-/im galaxy counts. Observed counts are taken 
from Stickel et al. (1998), Kawara et al. (1998), Puget 
et al. (1999) (FIRBACK), Oliver et al. (2000) (ELAIS), 
Juvela et al. (2000) and results of the fluctuation analy- 
sis of the Lockman Hole (Matsuhara et al. 2000). Most 
of these data suggest that the slope of the number counts 
is very steep, i.e. dlog A^/dlogS < — 2.5. This steep 
slope index strongly supports the rapid evolution toward 
z ~ 0.75. This suggested evolution is quite consistent 
with the requirements of the CIRB spectrum. 

We note that the Evolutions 1, 2, and 3 show signif- 
icant difference at 0.1 — 0.01 Jy in the FIR wavelength. 
Especially, only Evolution 3 can reproduce the constraint 
of Matsuhara et al. (2000). Thus a huge database, such 
as the data of ASTRO-F all-sky survey, can discriminate 
this difference and flx the evolutoinary history at this 
redshift regime (0<z<0.75). 

4.2.3. Submillimeter galaxy number counts 

The bottom-left and bottom-right panels of Figure 8 
show the 450- and 850-/im model number counts and 
observed counts obtained by SCUBA. It is now well ac- 
cepted that the sub-mm source counts are suitable for 
examine the status of the high-redshift {z > 1 — 2) 
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Fig. 8. — The multiband galaxy number count at infrared - sub-mm wavelengths. Same as in Figure 6, dotted lines 
represent the no-evolution prediction. Dot-dashed lines are the counts with Evolution 1, dashed lines are those with 
Evolution 2, and solid lines, with Evolution 3, respectively. 
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galaxies. The observed 450 fiia counts are recently ob- 
tained by Blain et al. (2000), and two upper limits are 
based on Small et al. (1997) and Barger, Cowie, Sanders 
(1999a). The symbols represent the observed 850 fun 
counts taken from Small, Ivison, Blain (1997), Hughes 
et al. (1998), Blain et al. (1999), Eales et al. (1999) and 
Holland et al. (1999), and the hatched region, Barger et 
al. (1999a). 

The 450-^m counts of Blain et al. (2000) can be re- 
produced by Evolution 3, because this evolutionary his- 
tory has a long plateau toward z = 5. On the other 
hand, the 850-/^m counts are somewhat problematic. If 
we attach importance to the FIRAS CIRB spectrum, our 
model slightly but significantly underestimates the source 
counts at 850 /im. We should consider the discrepancy 
found in our model predictions and detected sub-mm 
source counts, because in contrast, our model success- 
fully reproduced the FIRAS CIRB sub-mm spectrum. 
Since, in our estimate, IR galaxies are fewer at high red- 
shift compared with sub-mm-weighted model (e.g. Tan 
et al. 1999), this result seems qualitatively natural. 

Here we consider the difficulty in measuring a source 
flux from noisy data. Hogg, Turner (1998) clearly sum- 
marize and discuss the upward bias of the flux estima- 
tion in the study of source number counts. Very recently 
Hogg (2000) extensively examined the confusion effect on 
the astrometric and photometric measurements. Hogg 
showed an important result that if the observed num- 
ber counts has a steep slope like IR and sub-mm source 
counts, flux estimation must suffer from serious confusion 
noise, and most of the faintest sources can be spurious. 
According to Hogg's discussion, the sampling strategies 
adopted in the field of sub-mm survey are often insuf- 
ficient. Thus wc do not make attempts to search the 
solution which satisfy the present sub-mm source counts 
in this study. In order to settle down this problem, a 
much better rule-of-thumb for the same flux limit level 
for the confusion is required. Eales et al. (2000) also per- 
formed a Monte Carlo simulation of the source extraction 
from the sub-mm images, and clearly showed that the es- 
timated source fluxes are boosted upwards significantly, 
by the heavy confusion. 

We discuss another possible source of the sub-mm dis- 
crepancy. The hot gas in X-ray cluster of galaxies pro- 
duces a spectral distortion of the CMB: the Sunyaev- 
Zel'dovich (SZ) effect (Sunyaev, Zel'dovich 1972). At 
wavelengths shorter than 1.38 mm, the SZ effect appears 
as a positive source in the CMB and has a maximum 
intensity at 850 /im. The number count and redshift 
distribution of clusters are actively studied by many au- 
thors (e.g. Barbosa et al. 1996; Kitayama, Sasaki, Suto 
1998). Kitayama et al. (1998) prediced the sub-mm clus- 
ter number counts based on the Press-Schechter formal- 
ism (Press, Schechter 1974). The contribution of SZ clus- 
ters to the sub-mm numbers count is not negligible com- 



pared with sub-mm galaxy counts, while their contribu- 
tion to the CIRB is small (Kitayama et al. 1998). In ad- 
dition, SZ cluster is a diffuse source, and thus they might 
affect the flux estimation of sub-mm point sources in the 
cluster regions. Hughes, Gatztafiaga (2000) has already 
included the effect of the SZ clusters in their sub-mm sky 
image simulations. 

In addition, we should take the possible AGN contri- 
bution to the sub-mm source counts into account, though 
there are some observational constraints on their fraction 
in the sources (e.g. Haarsma, Partridge 1998). 

Instruments with a much better angular resolution in 
sub-mm wavelengths, like forthcoming Large Millimeter 
and Submillimeter Array (LMSA) and further, the in- 
ternational coUboration of large arrays, Atacama Large 
Millimeter Array (ALMA) will help to solve the confusion 
problem, and serve much information of the high- redshift 
status of galaxies. We stress that large- area surveys tar- 
geted at bright submillimeter sources are also necessary 
to fix the evolution at 1 < z < 2 (for the detailed exami- 
nation of this issue, see Takeuchi et al. 2000c and Hughes 
2000). 

5. Discussion 

5.1. Implication: High or Low Redshift? 

Our analysis shows that the steep slope of the IR num- 
ber counts and the peak intensity of the CIRB is at- 
tributed to relatively low-redshift (0.5 < z < 1) galaxies. 
Previous studies of the redshift distribution of IR galax- 
ies based on the CIRB spectrum and IR number counts 
tend to claim that the IR - sub-mm observations support 
a model of high star-formation rate (SFR) at large red- 
shift (z > 2), and consequently, imply a large number of 
heavily obscured distant galaxies. We focus on this issue 
and examine our results in this subsection. 

First we took care of the CIRB intensity at 140 /im 
(Hauser et al. 1998). As we discussed in Section 3, the ob- 
servation and calibration of COBE DIRBE at this wave- 
length is reliable enough. Therefore it seems strange that 
most of the previous studies completely ignored this ob- 
served point but focused on the other points. Moreover, 
some previous model predictions failed to reproduce the 
overall spectral shape of the CIRB. If we use normal dust 
temperature Tdust = 30 ^ 40 K, the spectral shape of the 
CIRB naturally requires a significant contribution of low- 
z IR galaxies. 

Gispert, Lagache, Puget (2000) coherently analysed 
the cosmic background from optical to millimeter wave- 
lengths and derived the cosmic SFH, and also found a 
rapid increase of cosmic luminosity density, and conse- 
quently, the cosmic SFR density. They also found that 
the contribution of the low-z IR galaxies to the CIRB 
around the peak is almost 100%. Their evolution of the 
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cosmic luminosity density at z < 1 is quite consistent 
with our result. Eales et al. (1999) studied their sub-mm 
survey and presented a similar SFH, although their peak 
of the SFH is located on slightly higher redshift. 

Recently Scott et al. (2000) performed a sub-mm 
follow-up observation of 10 known FIRBACK 170-/im 
sources and reported some detections at 450 and 850 ^m. 
They found that the sub-mm fluxes of the FIRBACK 
sources are fairly weak. We see, in their Figure 1, that 
the SEDs of these sources are consistent with those of IR 
galaxies at < z < 0.5. Scott et al. (2000) mentioned 
that there remains a degree of freedom to hypothesize 
the very hot Tdust- They carefully considered another 
redshift estimation, variation of dust emissivity and dust 
temperature, and concluded that these sources have red- 
shifts of < z < 1.5, which well agrees our result. 

Juvela et al. (2000) reported a result of galaxy count as 
a part of their ISOPHOT CIRB project. They performed 
multiband surveys at wavelengths 90 /zm, 150 /xm, and 
180 /zm, and relatively high source surface density sug- 
gested strong galaxy evolution. They compared the ob- 
served SEDs of the newly detected ISO sources with 
the IR-luminous galaxy SEDs of Lisenfeld, Isaak, Hills 
(2000), and found that the SEDs of the detected ojbects 
are consistent with those of dusy galaxies at redshift 
0-5<z<l. This also supports our estimated evolution 
of IR galaxies. 

We must note that, in addition, the density parame- 
ter of the Universe fio = 1 often assumed in many pre- 
vious models. Large flo suppress the faintest slope of 
the galaxy number count and intensity of the CIRB (e.g. 
Takeuchi et al. 1999). Therefore, a group of the models 
which prefer very high redshift IR galaxies can be rec- 
onciled with the limits of the FIRAS CIRB spectrum by 
the large cosmic density. On the other hand, the cosmo- 
logical constant Aq does not significantly affect the model 
predictions of the CIRB and IR galaxy number counts, 
compared with the effect of Hq. When we use the low- 
rio cosmology (0.2 < H. < 0.3), which is now regarded 
as highly plausible, the CIRB spectrum turns out to be 
a meaningful observational constraint of galaxies at high 
redshift. 

5.2. Cosmic Star Formation History 

In this subsection we discuss the star formation rate 
per comoving density as a function of redshift. This is 
generally called the cosmic star formation history (SFH). 
Since, as we mentioned in Section 1, the FIR luminosity 
of a galaxy is related to its star formation rate (SFR; e.g. 
Kennicutt 1998), we make attempts to convert the FIR 
luminosity density derived in the previous sections to the 
SFR density in the Universe. 

First, we estimate the comoving luminosity density 
in the FIR range of 40-120 fim, pfir(40-120). Using 
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Fig. 9. — The evolution of the far-infrared luminos- 
ity density within the wavelength range from 40 /im to 
120 /im, pfir(40-120). The integrated luminosity den- 
sity is derived from the monochromatic luminosity den- 
sity in two IRAS wavebands, peo/^m and piooAim with for- 
mula of Helou et al. (1988). 



the well-accepted formula presented in the Appendix of 
Helou et al. (1988), /9fir(40-120) is approximated as 



-19 



Pfir(40-120) = 3.26 x 10 

[2.58p^(60) -I- l.OOp^(lOO)] [Lq Mpc" 



(11) 



where Pi,{60) and Pi,(100) are the comoving luminos- 
ity density per unit frequency at 60 /xm and 100 fim 
[erg s~^Hz~^Mpc^'^], respectively. Helou et al. (1988) 
asserted that the relation in their Appendix is valid in 
the range of dust temperature Td and emissivity index 
7 expected for a galaxy. Thus, we expect that equation 
([n] ) is applicable to the FIR luminosity density of the 
Universe, since the FIR emission is originating from dust 
in galaxies. We show the history of pfir(40-120) in Fig- 
ure 9. Here we note that both Pi,(60) and p,^(100) at 
redshift z are calculated from the LF and SEDs shown 
in Section 2. 

Next, pfir(40-120) is converted to the SFR per unit 
comoving volume, psFR- Recently, Inoue, Hirashita, Ka- 
maya (2000) derived the conversion formula from FIR 
luminosity to the SFR. Their formula is applied to the 
relation between pfir(40"120) and psFR as 

1 .. 31 2.4 X 10-i°(l-r/) 
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Fig. 10. — The star formation history derived form cos- 
mic infrared background and IR number counts. The 
sohd line represents the case where the factor Cfir for 
the conversion from FIR luminosity to SFR is fixed at the 
value corresponding to the solar metallicity. The dotted 
line presents the case Cfir is 4.1 times larger because of 
the inefficient conversion from stellar luminosity to FIR 
luminosity. 



xpfir(40-120) [Lq Mpc 
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(12) 



where / is the fraction of ionizing photons absorbed by 
neutral hydrogen (i.e. (1 — /) is the fraction of ionizing 
photons absorbed by dust grains), e is the efficiency of 
dust absorption for nonionizing photons, and rj is the 
cirrus fraction of observed FIR luminosity. Since Inoue 
et al.'s expression is based on the luminosity in the range 
of 8-1000 /xm, we have converted it to the luminosity in 
40-120 Atm by the factor of 1.4 (e.g. Buat, Xu 1995). For 
convenience of discussions, we define the factor Cfir as 

Cfir [M, L- yr-] ^ M^il^:!!(i^ , (13) 
MR L o J J 0.4-0.2/ + 0.6e ' ^ ' 

or equivalently 

PSFR [M© yr-i Mpc-3] = 

Cfir • Pfir(40-120) [Lq Mpc-^] . (14) 

The above parameters / and e depend on dust-to- 
gas ratio, which is related to metallicity (e.g. Hirashita 
1999). Thus, Cfir depends on metallicity (Hirashita, 
Inoue, Kamaya 2000) and it varies as the metallicity 
evolution of galaxies. As the metallicity evolves from 



0.01 to 1 times the solar value, Cfir changes from 
5.8 X 10"^° Mq Lq^ yr-i to 1.4 x 10"^° Mq Lq^ yr"!. 
This is because the conversion efficiency from the stellar 
luminosity to the FIR luminosity increases as the dust- 
to-gas ratio (metallicity) becomes large. For simplicity, 
the cirrus fraction is set constant as r; = 0.5, since as 
stated in Hirashita et al. (2000), it is difficult to relate rj 
to the metallicity from our present knowledge. 

In Figure 10, we show the derived cosmic SFH. The 
solid line indicates the history determined by using con- 
stant Cfir, where the value for the solar metallicity 
is adopted (i.e., Cfir = 1-4 x lO"!" [Mq Lq^ yr"^])- 
The dotted line represents the history by using Cfir ~ 
1.4 X 10-1° 
5.8 X 10-1° 



[Mq Lq' yr- 



for 
for 



z < 1, while Cfir ~ 
[Mq Lq^ yr^i] for z > 1 to examine the 
extremely metal-poor case. In order to be consistent, we 
should model the cosmic chemical evolution. However, 
since the model produces further undetermined parame- 
ters, we do not treat such a model in this paper. We left 
detailed and consistent modeling to the future works (for 
reference on the cosmic chemical evolution, see e.g., Pei, 
Fall, Hauser 1999). We should also decrease rj as the FIR 
luminosity density increases, since the contribution from 
the cirrus component becomes smaller as the infrared lu- 
minosity becomes larger (e.g., Beichman, Helou 1991). 
If 7] is set as for an extreme case, Cfir becomes two 
times larger than that estimated above. Accordingly, the 
uncertainty of the Cfir induced by the cirrus fraction is 
estimated to be a factor of two. 

In a realistic situation, we expect that the cosmic SFH 
lies between these two lines in Figure 10. In summary, we 
conclude that psFR increases to ^ 10^^ Af© yr^^Mpc-'^ 
at z ^ 1 and remains almost constant or declines grad- 
ually for z > 1, very similar to the SFH presented in 
Figure 9 of Steidel et al. (1999) and that in Figure 6 of 
Gispert et al. (2000). 

6. Summary and Conclusion 

We estimated the galaxy evolution from infrared (IR) 
galaxy number count and cosmic infrared background 
(CIRB). We constructed an empirical galaxy count model 
based on the spectral energy distributions (SEDs) and 
luminosity function (LF), both of which are constructed 
from the IRAS results. The SEDs are made based on the 
IRAS color-luminosity relations and a tight far-IR-radio 
correlation. Using this model we estimated the plausible 
evolution of galaxies from the data in a nonparametric 
way. 

We found that a violent evolution of galaxies as a whole 
in relatively low redshifts (z< 0.5-0.75) should exist to 
explain the observed number count of galaxies in the IR 
wavelengths and the CIRB. Thus we conclude that the 
forthcoming galaxy surveys in the IR by future facilities 
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(e.g. ASTRO-F, SOFIA, SIRTF, FIRST, and NGST) 
have a crucial importance to understand the global his- 
tory of galaxies, since the IR wavelength is suitable for 
the investigations of the redshift range < z < 1. 

On the other hand, at present, star formation history 
at high redshift {z > 1) requires further observations 
and examinations of submillimeter (sub- mm) sources. In 
addition, source confusion severely affect the flux esti- 
mation in sub-mm. Instruments with much better angu- 
lar resolution in sub-mm wavelengths, like forthcoming 
LMSA and ALMA will help to solve the confusion prob- 
lem, and serve much information on the high-redshift sta- 
tus of galaxies. We stress that large-area surveys targeted 
at bright sub-mm sources are also necessary to fix the 
evolution at 2; > 1 — 2. 

Then we converted the FIR luminosity density to the 
SFR density in the Universe. We estimated the comov- 
ing luminosity density in the FIR range of 40-120 /im 
by the approximate formula of Helou et al. (1988), and 
Pfir(40-120) is converted to the star formation rate per 
unit comoving volume, psFR- We used the conversion for- 
mula from FIR luminosity to the SFR derived by Inouc 
et al. (2000), which took the following three parameters 
into account: the fraction of ionizing photons absorbed 
by neutral hydrogen /, the efficiency of dust absorption 
for nonionizing photons e, and the cirrus fraction of the 
observed far-infrared luminosity rj. Since / and e de- 
pend on dust-to-gas ratio which is related to the metal- 
licity, we can include the effect of the metal evolution in 
galaxies. We observe that the SFR density increases to 
^ 10^^ MQyr^^Mpc^'^ at z ^ 1 and it remains almost 
constant for z > 1, very similar to the star formation 
history presented by Steidel et al. (1999) and Gispert et 
al. (2000). 
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Appendix. Performance of the IRIS Far Infrared 
All-sky Survey 

In this Appendix we present the performance of the 
Japanese ASTRO-F (IRIS) all-sky survey. The bandpass 
system of the FIR instrument, ASTRO-F Far Infrared 
Surveyor (FIS) consists of two narrow bands, N60 (50 — 
70 ^m) and N170 (150 - 200 fim), and two wide bands, 
WIDE-S (50 - 110 Aim) and WIDE-L (110 - 200 ^im). 



We show the expected number counts at each band in 
Figures 11 and 12. The vertical dot-dot-dot-dashed line 
presents the 5(T-flux detection limit in each band. The 
detection limits are estimated as 39 niJy and 110 niJy 
for N60 and N170, and 16 mJy and 90 mJy for WIDE- 
S and WIDE-L, respectively (Takahashi et al. 2000). In 
order to obtain realistic galaxy counts, the wide bandpass 
wavelength range is involved in calculation. Figures 11 
and 12 show that several x 10^ galaxies per sr will be 
detected by the sets of narrow bandpasses, and roughly 
10® galaxies per sr will be detected in WIDE-S and -L. 

In order to determine the effect of galaxy evolu- 
tion from survey data, a significant sky area should be 
scanned, to suppress the variation of galaxy surface den- 
sity. Since galaxies are clustered on the sky, the nominal 
error bar estimated from the Poisson assumption is an 
underestimation, and we must take the galaxy angular 
correlation function into account. Consider a survey area 
fl, and divide it into small cells {Ailj} so that the num- 
ber of galaxies in the cell {Aili}, TVi = or 1. We set the 
mean galaxy surface number density J\f. Then we have a 
mean number in a solid-angle cell Af2i, A/i, 



{M,)=MArt, 



(Al) 



By definition, we have (Afj) = {Uf) = {Uf) = ■■■. We 
observe 

{KNj) ^ N^{i + w{e,j)) An.A^j (A2) 

((M-(M))(AG--(A/',))) 

==(A/-AG-)-(AA,)(AG-) 

= N^{1 + w{e,j)) AOjAOj - 7V^ AfljArjj 

^ M^w{e,j) AVl,AQ.j (A3) 

where w{9) is the angular two-point correlation function 
of galaxies. Next we consider the number of galaxies in 
the survey area O, N . We have 

{N) = Y^ (A/",) = I MdVL ^ JVn , (A4) 

(N^) ^ (J^AA.^AG) - 5](A/-2) + 5](A/-AG) 



i j i i^j 

i i^j 



AAdfi+ / / AA2(l + u;(6'i2))df7idfi2 



^N'n+Af^n^ +J^ / u;(6'i2)dr2idf72 . (A5) 
J Jn 

and thus 

{{N - (N))^) ^ (N^) - {Nf 

= Afn+Af^ f f w;(6ii2)df}idf72 ■ (A6) 
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Fig. 11. — The expected number counts in the ASTRO-F far-infrared all-sky galaxy survey at two narrow band 
filters, called N60 and N170. The wavelength ranges are 50 - 70 fim (N60) and 150 - 200 ^m (N170). 



If we assume that we take a sufficiently large area than 
the coherent scale of angular clustering, we can approxi- 
mate the above expression as 



{(N-{N)f)^Afn- 



In 



--Afn [i + Af w(e)dn 

Jn 



(A7) 



Here we evaluate the 'signal-to-noise ratio' of the number 
count S/N: 



S/N = 



(N) 






^j\fn{i + Afj^^w{9)dn) 



n 



J,,w{9)dn- 



(A8) 



We see that the S/N depends on the angular correla- 
tion strength and the solid angle of the survey, and is al- 
most independent of the surface density of the sources M. 
This fact shows that, in order to determine galaxy evo- 
lution from number counts, large-area survey is substan- 
tially required. The angular correlation function of the 
IRAS galaxies is w{0) = {9/0o)~°'^^ {Oo = 0?11: Lahav, 
Nemiroff, Piran 1990). The angular correlation of galax- 
ies detected by the ASTRO-F survey is obtained by using 



the scaling relation oiw{9) with detection limit, through 
relativistic Limber's equation (e.g. Peebles 1980). If we 
would like S/N = 10, a survey with ~ 50 - 300 deg is 
required, depending on the considered flux. We see that, 
although the deeper is desirable, but the wider is nec- 
essary. The ASTRO-F all-sky survey promises a unique 
opportunity to execute such a large-area survey. Thus we 
can expect for ASTRO-F survey to distinguish between 
Evolutions 1,2, and 3 clearly, and to fix the evolutionary 
history of IR galaxies up to z '-^ 1. 
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